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Abstract 

This note provides an explicit parametrization of all purifications of a mixed state in dimension 
2 and all joint purifications, if any, of two mixed states in the same dimension. The former 
is parametrized by SO(3, R), while the latter is parametrized by SO(2,R), except when the 
state being purified is already pure. Using this, we show how to calculate certain measures of 
quantum information and as a byproduct we show how to solve one variation of the classical 
Procustes problem. This manuscript was originally scheduled to appear on the arXiv on 
October 8th, 2002, but it did not due to alleged illegibility of the pdf version of the manuscript. 



1 Introduction 

The notion of purification of a mixed state plays an important role in several contexts, |IJ. 
It provides insight into the question of dccohcrcncc. It is important in quantum information 
theory from several points of view. For instance, many quantitative measures such as the 
maximal singlet fraction may be explicitly defined in terms of purifications. 

The purpose of this note is to provide an explicit parametrization of all possible purifica- 
tions of a mixed state in two dimensions, and joint purifications (if any) of two mixed states 
in a two dimensional Hubert space. In terms of this explicit parametrization, this note re- 
covers many of the quantitative measures mentioned above. In particular, it shows that the 
calculation of such measures reduces to optimization problems on the real orthogonal groups, 
SO(2,R) and SO(3,R). 

The balance of this note is organized as follows. In the next section a precise definition 
of what we mean by purification and joint purification is provided. The same section derives 
the parametrizations referred to before. The III section shows how to reduce calculations 
of measures in quantum information to optimization problems, and gives some instances of 
when this can be done in closed form. As a byproduct, some insight into solving certain 
variations of the classical Procustes problem is obtained. Section IV offers some conclusions. 
In an attempt to extend this mixed states of higher dimensions, we provide a Bloch sphere 
like characterization of three dimensional mixed and pure states in the appendix. The final 
section offers some conclusions. 

2 Parametrization of Purifications 

First, a mixed 2x2 state is psd, trace 1 matrix, p. A purification of p is a psd, trace 1, 
projection, P p , operating in C 2 <£> C 2 , such that the partial trace over the second C 2 factor 
of P p yields p. Likewise, given a pair, p\, p2 of mixed states in C 2 , every pure state P pitP2 in 
C 2 <g) C 2 whose partial trace over the first (resp. second) C 2 factor is pi (resp. p\) is said to 
be a joint purification of the pair p\, p2- 

Mixed states in C 2 may be represented in many fashions. However, one suitable choice 



for the purpose at hand is the following 
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i=l 

where /3, G i?, i = 1, . . . , 3 and i = 1, 2, 3 stands for i = x,y,z respectively. It is well known 
that || ||< 1, with equality precisely for those p which correspond to pure states. 

For the same reasons, a mixed state in C 2 <g> C 2 are best represented in the following 
fashion 

I 3 3 3 3 

P = ^ (-k + J! ^ CTi 72 + X 7l/2 CTi + H ^ CTfe ) ( 2>1 ) 
j=l t=l (=1 fe=l 

for real 0i,"fi,5ik- This representation is quite popular in the literature. However, to the 
best of our knowledge, it is never refined further to obtain a Bloch sphere like picture. This 
is needed for our purposes. The following characterization of mixed states and pure states 
follows from a direct calculation of squares of a Hermitian matrix (psd matrices are squares 
of Hermitian matrices and pure states are psd matrices which equal their squares) . To avoid 
circumlocution, we call the vector of Pi's as 0, the vector of ji's as 7 and the matrix of Ski's, 
5. 



Proposition 2.1 Every mixed state in C 2 ® C 2 is of the form in Equation (2.1), with = 
\{k0q + <5 7o), 7 == |(«7o + ^o A))i s = 5(^0 - (adj S ) T - 0qJq), for any /? ,7o € -R 3 , £>o € 
<?i(3, i?), K £ R satisfying || O || 2 + II 7o || 2 +Tt(5q'5o) < 4 and « is the positive square root 
of the difference of the RHS and the LHS of this inequality. Every pure state is of the form in 
Equation (Jy|), with = 5y, 7 = S T 0, S = -[adj (<5)] T +/3 7 T , \ \ || 2 + || 7 || 2 +Tr(S T S) = 3. 

Proof: The proof is a straightforward calculation. We will just record one important 
calculation going into the verification of this proposition, which will be needed for other 
purposes in this work. If pi,p2 are two mixed states represented via the form in Equation 



(2.1), then the trace of their product is 

Tr (p lP2 ) = <0u0 2 > + < 7i,72 > +Tr (S^S,)) (2.2) 

Next, note that the partial traces of such a p (pure or impure) are precisely the 2x2 
matrices \{^2 + 5Z»=i Pi a i)> 3(^2 + S*=i 7i <T i)- This, of course, directly implies that for any 
mixed state the lengths of and 7 is at most one. Further, a simple calculation, left to the 
reader, shows i) det(6) —\\ || 2 —1 for a pure state; ii) || \\ = \\ 7 || for a pure state. 



Returning to pure states, it follows from 5 = 
f3 || 2 -1, that 



-[adj (6)} T + /3 7 T , 7 = 8 T (3 and det(S) =\\ 



55 T = (l- \\ p \\ 2 )I 3 + f3(3 T (2.3) 

The following result says that this condition is essentially sufficient to determine purifica- 
tions of the state, \(Ji + 53-t=i 0i a i)- This result also provides a complete parametrization of 
such purifications. 

Theorem 2.1 Let \{h + 5Zi=i be a 2 x 2 mixed state. Then all possible purifications, 
P p may be parametrized as matrices of the form in Equation ( |2.l| ), with (3 the given the j3, 
7 = S T j3, S any solution of the system of equations: 

8S T = (l-H/3 | \ 2 )I 3 +(3f3 T (2.4) 
det(S) = ||/3 || 2 -1 

This system of equations is always solvable. Further, the general solution to this system (and 
thus the general purification, P p ) is provided by S = 5S, S € SO(3, R), with 6 one particular 



solution of this system, Equation (2.4). Further, the set of purifications of p is parametrized 
by SO(3, R) when || \\< 1 and by the unit sphere, S £ R 3 when || (3 \\= 1. 



Proof: First, it is clear that any purification has to satisfy the the system Equation (2.4). 
Proving the converse statement requires first proving that the respective system does have a 
solution, for any given (3 within the closed unit sphere in R 3 , and then that with the choice of (3, 
7 and S in the statement of the theorem, P p is indeed pure. In other words, this choice of (3, 7, 8 
indeed satisfies the defining relations for pure states, viz., || (3 || 2 + || 7 || 2 +Tr(8 T S) = 3, 
(3 = Sj; 7 = 5 T (3 and finally, 5 = -[adj (5)} T + /3-y T . 



Case I: || /3 ||= . In this case (3 is 0. So the system Equation (2.4) reduces to 

DD T = I 3 , det (D) = -1 

Clearly any matrix in 0(3, R) with determinant, —1, is a solution. In this case, 7 = D T (3 = 0, 
and Dj = DO = = j3. So of the defining relations for pure states, only the first and the 
fourth need checking. The first reduces to verifying Tr [DD T ] — 3 (since (3 = 7 = 0), which 



obviously holds. The final equation now becomes D = — (adj (D)) T , which also holds since 
det (D) = -1. 

If D G 0(3) is one solution to Equation ([O]), then so is DC, C £ 50(3, R). Conversely, if D 
is a second solution, then C = D~ 1 D exists and satisifies, 

CC T = D^ 1 DD T D~ T 1 = {DTD)- 1 = I 3 

Obviously det (C) = 1. So C e S*0(3, R). 

Case II: <|| (3 ||< 1: First, the matrix (1- || (3 || 2 )/ 3 + (3(3 T is positive definite. Indeed 
v T [(l- || (3 \\ 2 )I 3 +f3f3 T ]v = v T v-(3 T (3v T v + (< (3, v >) 2 is positive, for v ^ 0, since || (3 ||< 1. 
So there is always a matrix 8 satisfying 88 T = (1— || (3 \ \ 2 )Iz + f3(3 T ■ 
Let us compute the determinant of such a 6. We get 

(de^)) 2 = (l-H/3 || 2 ) 2 

Indeed, the eigenvectors of (1— || (3 || 2 )i3 + (3(3 T are (3 corresponding to eigenvalue 1 and 
any two vectors orthogonal to (3 (in i? 3 ) corresponding to the repeated eigenvalue 1— || (3 || 2 . 
From this the previous equation follows trivially. So, det(S), for a given solution S of the 
first equation in the system Equation ( |2.4| ), is either the positive or negative square root of 
(1— || P || 2 ) 2 - To ensure the negative square root, we multiply the given solution S by — 13 if 



needed. This also solves the first equation in Equation (2.4) and has the desired determinant. 

Defining 7 = 6 J (3, it is easy to verify (3 — 8^ and 8 = —(adj (8)) T + /3j T - Indeed, 
8j = 88 T (3 = (1— || (3 || 2 )/3 + (30 1 ' (3 = (3, yielding the desired conclusion. To verify, || [3 || 2 
+ II 7 II 2 +Tr((5<5 T ) = 3, we note first that upon taking trace on both sides of the first line in 
the system (|J) yields, Tr (88 T ) = 3 - 2 || (3 || 2 . Since Tr [8 T 8) = Tr (88 T ), it suffices to 
show that || S T f3 || 2 =|| (3 || 2 . But || 8 T f3 \\ 2 = f3 T 88 T (3 = /3 T [(1- || f3 \\ 2 )h + [3f3 T }f3 =|| (3 || 2 . 
Finally, to show 8 — — [adj (8)] T + (3j T , we will verify the transposed version. Since, 8 is 



invertible, we may premultiply both sides of the first line in the system Equation (|2.4|) by 



8 . This, bearing in mind the second line of the system (2.4), yields 



8 T = -[adj (8)} + 8~ 1 /3/3 T = -[adj (8)} + 1 p T 

since it was just shown that = f3 holds. 

To show that the 8 of all purifications is given by 8C, C E SO(3, R), with 8 one particular 



solution of the system (2.4), note first that 8 = 8C trivially satisfies (2.4). Conversely, if both 



(5 and 6 satisfy Equation (2.4), the polar decomposition theorem plus the fact that both S 
and S have the same determinant implies that there is a C £ SO(3, R) such that S = 8C (for 
an argument which eschews the polar decomposition theorem see the remark following the 
proof). 



Case III: || /? ||= 1 - In this case the system fl2.4|) reduces to 

DD T = P(3 T , det(D) = (2.5) 
Since per the first equation DD T is rank one, the second equation is superfluous. Once 



again there is at least one solution to Equation (2.5), viz 8 = [3(3 , for the given (3. 



Next, to verify that any solution to Equation (2.5), together with the given (3 yields a 
purification, we first observe that Sj = SS T (3 — (3, since || (3 || = 1. Just as in Case II, 
II 7 ll = ll ft ll = 1- This together with the obvious property that Tr (S T S) =|| (3 || 2 = 1 yields 
II ft II 2 + II 7 II 2 +Tr (S T S) = 3. To verify, the remaining condition, first note that the matrix 
S is also a rank one matrix. Indeed, the rank of 5 is the same as that of 88 T (this is valid for 
any square matrix). So S may be written in the form 5 = vu T for some vectors u,v G -R 3 . So, 

Tr (SS T ) = Tr (vu T uv T ) =|| u || 2 || v || 2 = Tr (/3/3 T ) =|| f3 || 2 = 1 

Hence, it holds that || u \\ || v \\= 1. So, dividing v,u by their lengths, if needed, it follows 
that v, u may be chosen to be of length one. Now, 

55 T = vv T = (3(3 T 

So either v = (3 or v = —(3. Abosrbing the negative sign if needed into u, we see S = j3u T for 
a length one vector u. Since u and f3 have length 1, and the group 50(3, R) acts transitively 
on the sphere in i? 3 , it follows that there is some C T £ SO(3,R) such that u = C T f3. So 
comparing S = (3(3 T with S = (3u T , we see i5 = SC. Since C £ S0(3,R), this verifies the 
claim. 

Finally, the assertion about the parametrization follows, since, when || j3 \\< 1, two purifi- 
cations with distinct S matrices are also distinct. When || [3 \ \= 1, however, two purifications 
are distinct only if C(3 ^ [3. This means the redundancy in the parametrization consists of 
the isotropy subgroup of SO(3, R)'s action at the point /3, which implies the stated condition 
on the parametrization in this situation. 



Remark 2.1 Alternative argument for part of Case II: In the following a different proof, 
which avoids the polar decomposition theorem, in Case II of the previous proof is given. This 
calculation may be of interest in its own right. Suppose S and S are two solutions of Equation 



(|2l|), then letting C = (S)~ 1 S, it is clear that det(C) = 1. Further, 

CC T = (S)- 1 ^- || p \\ 2 )h + A^K^)- 1 = (1- II P \\ 2 )[S T ~S}- 1 + where 7 = S T (3 

A direct calculation shows that S T 5 = (1— || f3 || 2 )i3 + 77 T - Writing the matrix on the RHS 
of this last equation as X + Y, we see that X and X + Y are both invertible and further Y 
is of rank one. By a trivial modification of the Sherman-Morrison- Woodbury formula, ||, it 
follows that (X + Y)- 1 = X- 1 - 1+Tr( 1 yx _ 1) ^" 1 y^" 1 , Applying this to X = (1- || (3 || 2 )/ 3 , 
Y = 77 T , yields 

ccf = ,1- 11 „ iraj^ft - T - 1 i fF ^] - ft 

i.e., C S 50(3, i?). 

Joint Purifications: Next, we suppose that two mixed states are given, i.e., pp — \{I2 + 
Ei=i ftcri), p 7 = i(/2 + X)i=i li a i)i with prespecified /3, 7 € i? 3 are given. When does there 
exist a pure state P in C 2 £g> C 2 , such that the partial trace of P over the second system yields 
pp, while that over the first yields p 1 . The aim is to parametrize all such Ps. 
Clearly, a necessary condition is that || (3 || = || 7 ||< 1. This is also sufficient. 

Theorem 2.2 Suppose /3, 7 € P 3 satisfy || /3 ||=|| 7 ||< 1. Then they can be jointly purified. 



Further, there is at least one solution, S, to the system Equation (2.4) which yields, per 
the prescription of Th 2.1, such a joint purification. Given one such solution, 5, the most 
general joint purification is given by SC, with C £ SO(3, R) satisfying C7 = 7. This is a set 
parametrized by SO(2, R), except when when (3 — 7 = 0, in which case it is parametrized by 
50(3, P) or when || f3 \\= 1, in which case there is a unique joint purification. 



Proof: Clearly, if we can find a solution 5 to the system Equation (2.4), which further 
satisfies the condition 5 T (3 = 7 for the given /3,7, the proof of Th 2.1 shows that the corre- 
sponding purification is indeed a joint purification. Suppose, for a specific solution, S sp , it 
holds that 6j p /3 7^ 7. Then, computing 



|| S sp [3 || 2 = (3 T 5 sp 5 T B = f3 T f3/3 T (3 = (3 T [l- \\ [3 \\ 2 }f3+ \\ [3 || 4 =|| [3 || 2 =|| 7 || 5 



The last equation follows from the hypothesis || (3 \ \ — \\ 7 ||. Thus, 5j p f3 has the same length 
as 7. Now, SO(3, R) acts transitively on spheres of any radius. So, there is a C T G SO(3, R) 
satisfying the condition C T 8^ p f3 — 7, i.e., 5 = 5 sp C provides one joint purification. Now of the 
purifications provided by S — 5C, C G SO (3, R), only those purifications which satisfy C T -f = 
7 will yield a joint purification. Further, via the same arguments in Theorem 1, the most 
general joint purification is necessarily supplied by C5 with C G SO(3, R) such that C7 = 7. 
The collection of all such C"s is, of course, the isotropy group at 7 of SO(3, R)'s action on this 
sphere. This isotropy group is conjugate to the isotropy at the vector || 7 || (1,0,0), which 
is precisely SO(2,R). Geometrically, all such C's are rotations in the plane perpendicular 
to the vector 7, while 7 is the axis of rotation. If || /3 ||< 1, then due to the invertibility 
of the S matrix of purifications, it follows that 50(2, R) paramterizes the collection of joint 
purifications. If || (3 ||= 1, then (3j T = [3{C^) T , for every C G SO(3,R) fixing 7. So there is 
just one joint purification. Finally, if 7 = 0, then the condition C7 — 7 is no constraint on 
CG SO(3,R). 

Remark 2.2 It is, of course, possible to induce on the set of purifications the additional 
structures in the orthogonal groups (or the sphere when || (3 || = 1). However, this may not 
be very useful. For instance, the Riemannian metric on the orthogonal groups may not be 
consistent with the any of the current notions of distance between pure states. However, in a 
certain sense, these additional structures will be employed later in this work. More precisely, 
in the next section, some calculations of quantum information measures will be reduced to 
optimization on the orthogonal groups. The fact that these problems have a solution follows 
from the compactness of these groups. Further, they can be reduced to optimization problems 
over products of closed intervals via Euler angles etc., 

3 Calculation of Certain Quantum Information Mea- 
sures 

With explicit parametrizations of purifications of a single mixed state and joint purifications 
of a pair of mixed states, it is possible to compute, either in closed form or via optimization 
over well defined quantities, several quantum information measures. Many such measures 



are often posed as optimization of scalar quantities over some pure states. Below we will 
give two examples where this can be done in closed form. The first is the maximal singlet 
fraction. A formula for this essentially appears in Q, where a full proof is not given. Further, 
the arguments involved in Q consist of reducing the S matrix of some mixed states into 
a normal form, which seems unmotivated. The proof below shows why that normal form 
naturally arises. Thus, the argument provided here may be seen as a complement to that in 
Secondly, we will compute the joint purification closest to a given ( impure) mixed state 
with the same partial traces. In general, this leads to a variation of the Procustes problem 
and this variation will be formulated as the solution of a concrete optimization problem over 
the interval [0, 2ir\. For the special case when this mixed state is the product state pp (g> p 1 it 
turns out that all joint purifications are at the same distance. We then explain this from the 
perspective of the functional being optimized in this generalized Procustes problem. 

Remark 3.1 While optimization over SO(3,R) or SO(2,R) may be viewed as constrained 
optimization problems and thus amenable to Lagrange mutiplier techniques, the methods 
used below avoid this. In part due to the nature of the function(al)s being optimized, it 
seems much better to use appropriate parametrizations of these groups and pass directly 
to an unconstrained optimization, than add further equations via the Lagrange multiplier 
method. 

Proposition 3.1 (see J^J) Consider a mixed state, p in C 2 <B> C 2 represented in the form 
given by Equation Denote the corresponding S matrix by 5 p . Then its maximal singlet 

fraction is given by |(1 + Yli=i a i)-> ^ det (S p ) < and by |(1 + a\ + a 2 — 03) if det (S p ) > 
Here, the Oi are the singular values of S pi with ui > <jj , i/ i < j, i, j = 1,2, 3. 



Proof: The maximal singlet fraction f(p) is defined via the equation f(p) = max < tp | p | 
ip >, where the maximization is over all pure states "0 which are maximally entangled. This 
collection of states is precisely the set of pure states locally equivalent to the Bell state. It is 
easy to see that this is precisely the collection of pure states, which when represented in form 



Equation <^1\) have ^ = 7 = 0, 5 = D,D e 0(3, i?), det (D) = -1. Now the following is true 

< i> I P I ip >= Tr (pp^) 



So finding f(p) amounts to maximizing Tr (pp^) over all which, when represented via 



Equation fl2.1]), verify (3 = 7 = 0, 6 = D, D G 0(3, E), det (D) = -1. 

Now by Equation ( |2.2| ), this is the same as maximizing the quantity j(l + Tr (<5jD), 
over Z? G 0(3, i?),dei (15) = — 1. Equivalently it is the maximization of Tr [(— 6j)V],V G 
50(3, R). This is closely related to the key step in the solution of the Procustes problem, 
g. Some of the steps require careful modification, since the solution in || uses optimization 
over the unitary group. In particular, the argument in (^] will not directly apply to the 
case det (5 P ) > 0. Therefore, only this situation is addressed here. If det (S p ) > 0, then 
det (-8J) < 0. Therefore, there exist S,T G 50(3, R) such that 

S(-5p)T = -diag (<j u cr 2 , cr 3 ) 

Hence we get, 

Tr [{Sj)V} - Tr [S T diag (a u a 2 , a 3 )T T V] 

which equals 

Tr [-diag (a 1 ,a 2 ,a 3 )T T VS T ] 

So this last quantity has to be maximized over V G 50(3,i?). Quite clearly the maximum 
occurs when T T VS T = diag (-1,-1,1). Now T, 5, diag (-1,-1,1) are all in SO(3,R), so 
such a V in 50(3, R) always exists and is unique. This then gives the stated expression for 
f(p). Note —I3 is not in 50(3,i?). So f(p) cannot be increased further. 

Next we look at the distance of the joint purifications of two mixed states given by Bloch 
vectors /3, 7 to a given impure density matrix, pp tl in C 2 <g) C 2 whose partial trace is also 
precisely the states represented by /?, a. The choice of distance is the Hilbert-Schmidt distance 
d(p,pp r/ ) 2 = Tr (p — Pf3 r( ) 2 - Suppose the 5 matrix of of P/3,7) is denoted E. Then from 
Equation fl2.2| ) it follows 

d(p, m7 ) 2 = i(l + Tr (6-E)(6-Ef) 

Fixing one choice, D for S it follows from Th 2, that the most general such S is given by DC, 
with C in 50(3, R) satisfying O7 = 7. So the above quantity becomes Tr(DD T + EE T — 
DCE T -EC T D T ). Since D, E are fixed minimizing this quantity is the same as maximizing 
Tr{DCE T + EC T D T ). But, EC T D T = (DCE T ) T . Hence, using the cyclic invariance of 



trace, the problem reduces to maximizing, over all C G 5*0(3, R) satisfying O7 = 7, the 
function 

Fd,e(C) = Tr (E T DC) (3.6) 
We begin with a simple observation 

Lemma 3.1 Suppose pp n is the tensor product of the mixed states corresponding to (3 and 
7, then the function Fd,e(C) is constant and equals || (3 || 2 , i.e., every joint purification is 
equidistant from /?/3 j7 . 

Proof: Since E = /3 7 T , in this case, it follows F DiB (C) = Tr {E T DC) = Tr (DCE T ) = 
Tr (D^f3 T ), since C7 = 7. Further, — (3 for pure states. So this reduces to Tr (/3/3 T ) =|| 
/?ll 2 - 

We now address the general situation. From this analysis a geometric interpretation for 
the previous lemma will emerge. Now maximizing Fd,e(C) would reduce to the key step 
in the usual Procustes problem, but for the restriction that C7 = 7. For 7 = this is no 
restriction, though as in the calculation of the maximal singlet fraction care has to be taken 
since the optimization is over 5*0(3, R). Therefore, we will study only the 7^0 situation. 

Proposition 3.2 The maximum of Fd,e{C) is given by the maximum of a differentiable 
function f(9) (defined in Equation fl3.7| ) below) over the interval [0, 2tt]. Thus, this maximum 
exists and can be found by comparing the values of / at 8 = Q and at the critical points of / 
in (0,2tt). 

Proof: By the singular value decomposition 

E T D = VT,W T 

where V is a real matrix with columns Vi eigenvectors of E T D(E T D) T , W is a real ma- 
trix whose columns, Wi are the eigenvectors of (E T D) T E T D and S is a diagonal matrix 
diag(ai, <T2, 03) where the Oi are the singular values of E T D. 
Once again by the cyclic invariance of Tr it follows 

F d ,e{C) = Tr (Y,W T CV) 

To find the diagonal entries of W T CV, we expand v i} Wi in some orthogonal basis, whose 
first member is member is 7 Denoting the components of Vi, Wi in the first direction by pi, Si 



respectively, we find 



Cvi = pij + R g xi,i = 



1,...,3 



Here pi =< v i: j >, while the Xj are precisely the orthogonal projections of the Vi onto the 
plane perpendicular to 7 (and thus,Xi is uniquely determined by the Vi). In fact, since 7 
is orthogonal to this plane, = Vi — Pij. Rg is the rotation through that the matrix C 
performs. 

Denoting by j/j, i = 1, . . . , 3 the orthogonal projections of the Wi onto the plane perpen- 
dicular to 7 and Si —< Wi, 7 >, we find 



So maximizing this function will yield the distance of pp i<7 from the set of joint purifications 
of the mixed states represented by the Bloch vectors /?, 7. 

Remark 3.2 The situation covered by Lemma can now be explained as follows. E — fi^f T ■ 
So the singular values of E T D are the eigenvalues of E T DD T E. This is the rank one matrix 
II 7 II 2 77 T - Hence the singular values are (|| 7 || 2 ,0,0). Further, v\ = w\ = 7. Hence, 
Xi = 0. So Fd^e{C) us independent of 0, i.e., of C. Finally, since || 7 || = || (3 ||, the stated 
value for F DtE {C) is indeed obtained. 

Remark 3.3 In Proposition 3.1, the optimization over SO(3,R) did not require more than 
elementary aspects of the group structure of SO(3, R). For related optimization problems we 
found the following parametrization of 50(3, R) useful. Let the normalized eigenvector of the 
generic D e SO(3,R), belonging to the eigenvalue 1, be written in spherical coordinates in 
the form e_o = (cos 0d sin 4*d , sin 0d sin 4>d , cos 4>d)- Denote the angle of rotation D performs 
(counterclockwise) in the plane orthogonal to this eigenvector be tpD- Then the effect of 
D on any vector in v e R 3 is < v, en > en + R^, D (v— < v,eo > £d)- Since the typical 
optimization problem arising in contexts similar to the ones in this work involve inner products 
and norms, this representation of 5*0(3, R) seems superior to others (over Euler angles or 
Givens rotations, for instance). In this representation the optimization reduces to optimizing 
a function F(</>,6,ip) over [0,tt] x [0, 2tt] x [0,2tt]. 



3 




(3.7) 



4 Conclusions 

This note yielded a complete parametrization of purifications and joint purifications of 2 x 2 
density matrices. This enabled a reformulation of the calculation of quantum information 
measures as optimization problems over the real orthogonal groups. In particular, a solution 
to one variation of the classical Procustes problem was provided. It would be interesting 
to extend this to density matrices in higher dimensions. The first ingredient in this a clear 
description of density matrices and pure states, going beyond the fact that they are expressible 
as real linear combination of certain matrices and necessarily have trace 1. In other words, a 
full characterization of this real vector of coefficients in this expression is desirable. This is 
partially addressed by the appendix. 
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6 Appendix 

How may one generalize this explicit parametrization of purifications? The first ingredient 
is a Bloch sphere like picture of density matrices in the appropriate C™. For this there 
are two standard points of departure beyond the setting here. One is to consider higher 
tensor products of two dimensional spaces or to consider twofold tensor products of spaces of 
dimension higher than two. In an attempt to achieve this for the second route, we consider 
the question of describing, in a Bloch sphere fashion, density matrices of a single system first. 
This already is quite an arduous job as may be seen below. 

First we represent a typical n x n density matrix in the following form: 

^ n 2 -l 

p =-(/*+ y\ m (6.8) 

n •f— i 

2—1 

where the matrices A^ satsify i) {i\k, k — 1, . . . , n 2 — 1} is an orthogonal basis for su(n); ii) 

n 2 1 

their "Jordan" commutator satisfies, (A&A; + A/Afc) = ^7„ + J2i=i dkuM, with the dku a 
symmetric tensor. Such bases always exist. In principle, the dun can be found, though we are 
aware of their enumeration only for n = 3, 4 for one particular choice of such a basis, S. For 



n — 3 this basis is precisely the set of the Gell-Mann matrices. The reasons for choosing this 
representation is twofold i) eventually we wish to take partial traces of density matrices in 
C n ® C n . Therefore, having a basis for density matrices in C n which is maximally traceless 
will facilitate this computation, ii) The A^'s properties are well known. 

For any n, using the symmetric tensor, ctyfc, define a new vector iUj £ R n ~ starting 
with two vectors x, y G R n _1 via 

n 2 -l n 2 -l 

x U y = ( ^2 d ljk Xjy k , . . . , ^ d ljk Xjy k , . . .) 
j,fe=i j,fe=i 

The goal now is to characterize the form of the vector (3 £ R n ~ which ensures p is a 

density matrix. Since p is the square of a Hermitian matrix, M, we take a generic such M, 

square it and insist its trace be 1. This leads to the following characterization of mixed states 

and in particular, pure states. 



Proposition 6.1 Every density matrix can be represented in the form in Equation (6.8) with 

(i = it A) + where & is an y vector in i? " 2 ~ 1 with II A) ll 2 < 7J r and K = +\/ " 2 ~ 2 l IA) " 2 ~- 

Conversely any Hermitian matrix admitting such a representation is necessarily a density 



matrix, p is pure precisely if it can be represented in the form in Equation 3.8) with < 

(3, (3 >= and (n -2)f3 = /3Uf3. 

Note that when n = 2, this is precisely the usual Bloch sphere, since the Pauli matrices 
anti-commute, i.e. (3 U (3 — if n is 2. For n — 3, states are pure if the vector (3 S i? 8 is of 
Euclidean length V3 and (3 = (3 U /3. One can now use this to write down a Bloch sphere like 
picture for density matrices in C n ® C™ . This characterization involves both the Lie product 
and the Jordan product. However, since even for single systems, the picture provided above 
requires further analysis the details will be considered in a future study. Indeed even analyzing 
the implications of an equation of the type (3 = /3U/3 for density matrices in C 3 requires further 
work. However, given specific Hermitian operators (for one or two particles) one can use this 
characterization to check if they are indeed density matrices and even if they are pure. 
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